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Abstract. We propose a geometrical approach to the investigation of Hamil- 
tonian systems on (Pseudo) Riemannian manifolds. A new geometrical criterion 
of instability and chaos is proposed. This approach is more generic than well 
known reduction to the geodesic flow. It is applicable to various astrophysical 
and cosmological problems. 

1 Introduction 

Our main aim is to study global behaviour of dynamical systems on (Pseudo) 
Riemannian manifolds using their local, geometric properties. The guiding prin- 
ciple comes from the well known example of dynamical system: geodesic flow. 

If geodesic flow on Riemannian manifold is given, then to work out principal 
features of the flow, so-called sectional curvature should be investigated. And if 
it turns out that sectional curvature is negative definite, then one deduces that 
geodesic flow is Anosov with hyperbolic geodesies. 

As is well known, any Hamiltonian system (with the Hamiltonian of the 
form kinetic plus potential) is reducible to a geodesic flow. Therefore, above 
mentioned procedure may be used for the new system to get some characteristic 
features of the original system. In our opinion it is not a proper way to study 
dynamical systems. First of all, it does not work for all Hamiltonian systems 
(there must not be turning points i.e. kinetic energy must be positive), then in 
order to reduce to a geodesic flow one has to change the metric, Levi-Civita con- 
nection, and time parameter. And since sectional curvature essentially depends 
on the above mentioned values, one must make inverse transforms, in order to 
obtain characteristic parameters for the original system, which is not a trivial 
problem for typical systems. 

Here, in particular, we are looking for a local, geometrical object (it should be 
the sectional curvature in the case of the geodesic flow) for Hamiltonian systems 
which should characterise instability, hyperbolicity of dynamical systems. 

We adopt the following approach to the problem, which may be considered 
as an alternative to the approach mentioned above. Instead of changing met- 
ric, Levi-Civita connection (covariant derivative), and time, we change only the 
derivative. We deflne a new covariant derivative -D„ for any flow, which is deter- 
mined uniquely and for simple cases, coincides with the well known derivatives 
(Vu-covariant derivative, .Fu-Fermi derivative). By means of the derivative we 
introduce a coordinate independent equation for the flow invariant separation 
vector (perpendicular to the velocity). For the geodesic flow it is nothing but 
the Jacobi equation. As a by-product we represent any invariant vector fleld as a 
sum of vectors, which are perpendicular and parallel to the velocity respectively. 

t Published in "Proceedings of the Fourth Monash General Relativity Workshop" Eds 
A.Lun, L.Brewin, E.Chow (Department of Mathematics, Monash University, 1993), p. 38 



Then we introduce a geometrical object flu negativity of which means hyper- 
bohcity of the corresponding flow. This method does work even for systems 
with turning points and moreover, which is crucial for the Einstein dynamics, 
the geometrical approach applies for Pseudo-Riemannian manifolds as well. 

2 Geometry of Dynamical Systems 
2.1 Dynamical Systems 

Let S' be a smooth vector field on the tangent bundle r : TM ^ A/ of a 
Riemannian manifold endowed with a Riemannian metric g, its associated Levi- 
Civita connection K with covariant metric V (cf. ^) 

S : TM TTM -.u^iu, -dV) G Hor{u) ® Ver{u) , 

where V : M ^ and dV is the gradient vector of V, i.e. vector field such 
that 



Hor(u) and Ver{u) are horizontal and vertical subspaces of T„TAf respectively. 
We denote by u{t) — /*uo an integral curve of the vector field S passing through 
the initial point mq G TM i.e. 

ii = S{u) , 

and m(0) = uq . The group of diffeomorphisms {/*}, t ^ R 

f : TM TM : uq ^ u{t) , 

is called a dynamical system or flow of the vector fleld S on TM. The curve 
u{t) — /*uo is called the flow line or integral curve starting at uq, c = r o u is a 
trajectory. 

Lemma 2.1 If u(t) = f^uo is an integral curve of the vector field S starting at 
Uq, then c{t) —TO u{t) is a solution of the following equation: 



determined by c(0) = t o uq and c(0) = uq. 

Lemma 2.2 // VcC = -dV, then c{t) = /*c(0). 

We are interested in the behaviour of the nearby integral curves of the in- 
tegral curve u{t) = /'mq. To deal with that problem, next we introduce an 
equation for the invariant vector fleld along the curve c{t). 

Lemma 2.3 Invariant vector fields Z{t) along the integral curve u{t) — f^ua, 
i.e. 



{dV,X) = g{dV,X) , 



yx e TM , 




(1) 



Z{t) = TfZo , 



Zo e Tu,TM 
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are m 1 : 1 correspondence with the solutions of the following equation: 

Vlz + ^uiz)~Hv{z)^0 , (2) 



where ^uiz) = Riem{z,u)u, Hv{z) — The correspondence is given by 

z{t) ^ Z{t) = (z(t), V„z(t)) e Hor{u{t))®Ver{u{t)) . 

Note that one might consider Eq. (0) as an equation obtained from Eq. |^ 
merely by covariant differentiation with respect to z. 

Thus we arrive at the foUowing system of variational equations 

c = u , 

Vuu = -dV , (3) 
V2z + 5R„(z)-i/y(z) - . 

It is worth mentioning that the vector z describes the motion of points (but not 
integral curves) of the nearby integral curves of u{t). However, in most cases it 
is important to know the behaviour of nearby integral curves (which means z is 
orthogonal to velocity u) with the same integrals. 

If = then the flow corresponding to S is called a geodesic flow and 
Eq. reads as 

c — u, 
V„u = , 
V2z + 3fi„(z) = , 

the last equation is the Jacobi equation. 

We say that the Jacobi field z (solution of the Jacobi equation) is orthogonal 
to u if the following conditions hold g(w, z) = and z{E) = g{u, V„z) — along 
integral curve. One can represent any Jacobi field as follows: 

z = n + [n{E)t + g{u, z)\q]u , (4) 

where n is orthogonal to u, i.e. g{u, n) — g{u, V„n) = . One may prove (see 
12) that in the case of the geodesic flow n is orthogonal to u for all t if it is so 
at t = 0. So if z{0) = n(0) and n{E) = then z{t) = n(t) for arbitrary t. It is 
not true for a general flow. 

It is an important property of the geodesic flow, because it means that the 
equation for n is the same as for z and one should only carefully choose the 
initial conditions. 

As is well known (see |3]) by changing the time parameter i, metric g, and 
linear connection K (covariant derivative V) one can reduce a Hamiltonian 
system to a geodesic one. 

Theorem 2.1 On TM\e = {m G TM \ g{u, u) = 2{E - V{t o u)) ^ 0} , there 
exist 

Gl. s = ip{t) parameter such that ip ~ \f2\E — V\ , 
G2. Metric g = \E - V\g , 
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G3. Linear Levi-Civita connection K with a covariant derivative 
VxY = VxY + Xi(b)Y + Y{^)X - g{X, Y)d4> , 

where (f> = ^hi\E — V\ , such that instead of Eq. 0) one has 

{i ^v , 
Vln + ^^{n) = , 

where 7' = Riem is the Riemannian tensor of the linear connection K , 

y{u,u)^E-V ^ g{v,v)^sgn{E~V) . 

So any Hamiltonian system can be reduced to a geodesic one. But this procedure 
has several demerits: 

1) time parameter t changes to afhne parameter s, to recover t one has to 
find function ip; 

2) metric g changes to g it changes Levi-Civita Hnear connection and the 
measure of length of vectors. 

But as is well known local stability/instability does depend on metric and 
time. And in order to reduce a Hamiltonian system to a geodesic one we 
have changed the both concepts. Therefore, to make a conclusion about stabil- 
ity/instability of original system one has to make the inverse transform. Which 
turns out to be a complicated one. To overcome the problem we suggest an al- 
ternative method. For this purpose we will introduce a new covariant derivative 
for dynamical systems and derive desired equation by help of the derivative. 

Now we will see how to overcome these difficulties. 

2.2 Covariant Derivative of Dynamical System 

Let us introduce a new covariant derivative along curve c. Given a smooth 
curve c: R^ M, ce C°°{R, M), with the velocity vector w = c = Tc. . We 
denote the set of all smooth vector fields (functions) along curve c by Xc {Tc) 
respectively. First we define the set of curves C along which one may define the 
covariant derivative. 

Definition 2.1 We say that c Cz C if there exists unique smooth (1, 1) type 
tensor field Qc along c such that 

Qc = — for all u such that g{u, it) 7^ . 

g[u,u) 

Clearly, if g(u, u) ^ along a curve c then c G C . Another important subset of 
C is established by the following Lemma. 
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Lemma 2.4 If c is a non-constant solution of the Eq. then c € C and 



Let Vc = S-Q,. 

Definition 2.2 Given c E C and X E Xc.. We will say that X E Xjr' i.e. X is 
a D-differentiahle along c if there exists unique smooth vector field D^X E Xc 
such that if g{u,u) ^ then 



Obviously, if g{u,u) ^ along c then X^ = Xc . If g(X,u) — fg(u,u), where 
f E Tc then X E X^ . In particular, u E X^ . Besides, if g{X,u) — then 
X E X^ . Du -.X^^Xc : X^ DuX is a derivative (see 0). 
Some easy-to-check properties of : 

1. DuX = QcCuX if QcX = X [VcX = 0). 

2. DuX = Vc^uX if VcX - X [QcX - 0). 

3. Du = Vu if VuU ~ i.e. c is a geodesic. 

4. Du = Fu if .g(u, u) = const ^ 0. Fu is the Fermi derivative. 

5. DuX = VcVuVcX + CuQcX . 

The derivative along c can be naturally extended from vector fields to arbi- 
trary tensor fields as follows: 

Dl. Du is a linear mapping of tensor fields of any type along c to tensor fields 
of the same type, which commutes with contractions; 

D2. Du{Ti ® T2) = DuTi ®T2 + Ti® DuT2; 
D3. Duf = u{f) = for any function /. 

Theorem 2.2 There exists exactly one derivative determined by |3P with 1-2 
and D1-D3 properties. 

Definition 2.3 Given c E C and X E Xc. We say that X is a D-parallel vector 
field along c if DuX — 0. 

Theorem 2.3 For any Y E Tc{a)M (if g{u, u)|o = then QcY — 0) there exists 
unique parallel vector field X along c such that DuX = and X{0) — Y . 




DuX = Vc^uVcX + QcCuQcX . 



(5) 
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If QcX\Q = = QcY\o and D^X = = D„y then g{X, u) = 0, g{Y, u) = 0, and 
g{X,Y) = const along curve. Notice that I?„u = 0. 

Now we derive the desired equation for n — VcZ by terms of the covariant 
derivative. 

Theorem 2.4 Given c which satisfies the Eq. 0]). Then any solution of the 
Eq. |0j has the following form (we consider limit-values of all objects atg{u,u) = 

-2n{V)+n{E) 



2{E - V) 

where n — VcZ satisfies the following equation 

n{E) 



-dt 



u , 



(6) 



Qcn\to 



E-V 
QcDun\t, = 



VcdV , 



where 



-Hv{n) 



2{E-V) 



dV 



E = ^g{u, u) + V — const along integral curve, 

n{E) ~ z{E) = g{u, V„z) + z(V) = const along integral curve. 

Let us mention that if ^ = 0, and E — ^ afong the integral curves (geodesies 
with g(u, u) — 1) then decomposition of z Eq. © coincides with the Jacobi field 
decomposition Eq. If we consider all integral curves of fixed energy then 
the following claim holds. 

Corollary 2.1 // Qcz\to = 0, z{E) = 0, then 



n{V) 
E-V 



dr 



where n = VcZ satisfies the following equation 

Dln + n^{n) = , 
Qcn\to = QcDun\to = . 

Corollary 2.2 If g{u,u) ^ (there is no turning point) and n{E) = then 



(7) 



Dln + nu{n) = Q , 
5r(u,n)|o = g{u,Dun)\(i = 



Vln + 'Ry{n) = , 
g(u,n)|o = V„n)|o = 



One can derive an equation for the length of n from Eq. Let n — where 
= £ and therefore g{^,£,) = 1 and g{S,,u) = 0. Then 



(8) 
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where 



K{(„ u) = g(5R„(0, C) Riem{u, ^, u, £,) , 



2{E-V) ' 



And 



Dl^ + 2XDuC + (A + A2)C + n^iO = , 

where A = £/£. 

It is worth noting that and Wu(C, ^) are smooth functions with respect 

to their aU arguments (a;,M, ^) along integral curves, even at turning points. 
And if u(0) = then 

Luoi^,o = i^^^uito = \\dv\m,o , 

where k is the second fundamental form of the V = E submanifold in M and 

{dv,o^o. 

Theorem 2.5 // SMe = {u e TM\^g{u,u) + V{t o u) = E} is a compact 
manifold, g{u,u) + ||rfV^|p ^ (i.e. there is no singular point), and 



(9) 



where u G SMe, g{S,,C) = I7 o,nd g{u,^) = then the dynamical system is 
Anosov. 

Condition © may be fulfilled unless dynamical system has an integral on SMe- 
One may average oJui^, by replacing ^ ^ with 



to get 

where 
and 



B{u) 



1 



d-i V 2{E - y ' 

hv{u,u) 3\uAdV\^ 



Ric{u, u) - + 



g{u,u) g{u,uY 



\X A y p - g{X, X)g{Y, Y) - g{X, Yf 



Bq(x) = lim B(u) 



d-1 



-AV- 



hv{dv,dv) 

\\dV\\^ 



One may speculate that a solution of the Eq. © in average is greater than a 
solution of the following "evolution equation" (cf. 0,0) 

e{t) + B{t)i{t) ^ , (10) 

where B{t) = B{u{t)). Notice that Eq. is an exact if d = 2. 
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2.3 Pseudo-Riemannian Manifolds 



The impossibility of direct application of the results developed for Riemannian 
manifolds is evidently connected with the Pseudo-Riemannian signature of the 
metric, e.g. it is easy to see that two close geodesies can diverge in Pseudo- 
Riemannian manifold, while the length of the separation vector may remain 
close to zero or even negative. Therefore not only new criteria are needed here 
but one should also redefine the stability properties themselves. 

Consider a dynamical system on Pseudo-Riemannian manifold M (see 
Let c be a trajectory passing through a point c(0) = m G M and {E^} a = 
1, . . . ,d basis for T^M. One can propagate {E^} along c so that 

DuEa = 0, EM^El, Va=l,...,d, 

and get a basis {Eg} along the curve c for any T^i^f-^M . We will call it I?-basis. 
Any vector X G Tc{t)M can be presented by means of the -D-basis {Eg} 

X{t)=X'^{t)Ea . 

The expression 

gEiX,Y)=Y,X''Y- 

a 

defines i?-metric on TM. Length of the vector X with respect to the basis {Ea} 
is defined to be 

\ a 

Let {Ea'} be another Z3-basis along the same curve c. Then a non-singular 
matrix $ exists, such that 

E^a=Y.^a''E^„. 

v 

So far as both {Ea} and {Ea'} being D-bases are D-parallel transported along 
c, this relation must hold for constant <!>. Therefore 

Ea=Y,^<^''' Eb' ■ 

v 

Thus we have 

a 

Since $ is non-singular then 

Y^^x^'f > , 

6' 




8 



forif =0 weget 

b' 



it is a contradiction with non-singularity of Therefore we have a positive 
defined matrix 

6' b' 

Space of aU X with ~ 1 is a compact space, therefore there are positive 

constants a and /? (depending only on basis) such that 



0<«||X|||< ||X|||, </?||X||| . 
For the vector field n: n — n°-Ea, DuU — h'^Ea, and 

eE{n)^{\\n\\l + \\D^n\\l)^/' . 

Now we can define stability, hyperbolicity, etc. of integral curves. 

Definition 2.4 Integral curve u is called to be linear stable if We > 0, 35(e) > 
such that £B(ri(0)) < 5{e) lEinit)) < e , for all t > 0, where n is a solution 
of the Eq. Q). 

// u is not linear stable it is called linear unstable. 

One can readily verify that the definition does not depend on the choice of 
D-basis {Ea}. Lyapunov characteristic exponents x for integral curve u is 
defined as follows 

/ N r ln^£;(n) 
X(u, n) = iim sup . 

t^oo t 

Evidently x does not depend on choice of {-Eq} either. One can give a definition 
of hyperbohcity as well (see |5]). 

Let n = where \\n\\E — t and therefore \\S]\e = 1, then (cf. Eq. ^) 

i+ [gEi^uiO.O - gE{DuLD^m = , 

where 

ig{i.dV)9E{LdV) 



dEi^uiO^O = 9Ei^u{0,0 - 9e{Hv{0,0 



2{E-V) 
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3 Einstein Dynamics on Superspace 

3.1 World, Universe, Superspace 

The set of d-dimensional Universes will be described as follows (see We 
assume, that the Universe is closed (compact and without boundary). 

By A4'^~^^ we denote the set of all d+ 1-dimensional, smooth (from the class 
C, r > 2), oriented, compact manifolds {d > 1), 

TVJ'^+i = {M'^+i} = { all d + 1-dimensional smooth, 

oriented, compact manifolds } . 

World (i.e. spacetime with material fields) will mean the following triad: 

iM^+\giM),^M)) , 

where M'^'^'^ G A4'^~^^, and g{M) is a smooth Riemannian metric on M''-'^^, 
$(M) is a smooth scalar field. Denote the set of worlds by W'^'^^, 

W^+' = M = {{M''+\g{M),^M))} = {(Af„,g^,$^)} . 

Let Tc{M''''^^) be a set of smooth functions on M'^'^^ without singular critical 
points (Morse's function) such that / e Tc{M'^+^) if 

/:M'^+i^5i = [0,2^]/{0,2^}, 
/ e C"'(M''+i) , 

if dM'^+^ ^ then f{dM'^+^) = d[f{M'^+^)] , 

where dN is the boundary of the manifold N . For every c £ , f £ J-'c(A/''+^) 
we denote 

MM''+']^{x\xeM^+\f{x)=c} , 

y{/,[M'^+i]} = {x\x e /,[M'*+i],d/|^^[M^+i](a:) = 0} . 

For given w G W^^^, f G TdM^), and c S S"^ we have the following triad: 

u{w,f,c) = (/c[M^],5,0) , 

where g is the metric induced on /c[M^,] hy g^, <j> = '^w\f^[M^]- Universes (i.e. 
space with material fields) are members of the set U''' , 

U U U"(^'/'^)' 

U''^{u}^{{Tu,guAu)} ■ 
According to Morse's theory 

^ T,'^ u n'^ U {0} , 
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where S'' is the set of all d-dimcnsional smooth, oriented, closed manifolds with a 
smooth Riemannian metric and a smooth scalar field on them. We will consider 
the empty set as a trivial manifold. 

At Morse's reconstructions fl'^ are critical sets with a given metric and scalar 
fields. 

In order to construct the set as a topological space one needs to introduce 
a topology, a system of open sets. We will do it as follows: first we define the 
set of "smooth" curves on W^. 

A mapping A : (0, 1) is called a "smooth" one in ZY<^, if 3w e 

and 3/ e Tc{Mw) such that 

X(^j){T) = u{w,f,T), Te(o,i)c5i. 

O is the strongest topology on the set W^, such that every "smooth" curves from 
the following set 

U U 

is continuous on (0, 1). 

By Suporspaco wc mean an endowed with the topology O {W^, O). We 
denote that space by again. 

The space of Universes with given manifold M we denote by Uf^, 

Uti = {uGU''\Tu = M} . 
3.2 Geometry of Superspace 

Let us fix any M G T,'^ and consider Uf^. In this case is the space of all 
smooth Riemannian metrics and smooth scalar fields on M. It is known that 
there exists smooth Banach structure on such spaces. 

We will introduce a metric on Uf^ such that kinematical part of the Hamil- 
tonian given by ADM formalism could be expressed by that metric. So we have 
the following metric on Uf^ 

g\g, </)](fc, x; h,6)= [ [-trik)tr{h) + tr{k x h) + x6]N-^dix{g) , 

where N is a positive function N : M ^ — {x e R\x > 0} and 

dii{g) = (det gf/'^dx'^ A . . . A dx'^ , 

tr{k) = g^^kab , {k x h)ab = kac9'"^hdb ■ 

The metric has an inverse metric and 
Q-^[g,4>]{'JT,p;Q,q) = (^-^^tr{'K')tr{Q') + tr{'K' x q') +pv) Nd^l{g) , 
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where 

(7r,p), (£», q) e T^g^^)l^M 1 = 7r'd^(.g), ...,q = q'dn{g) . 

By means of this metric one can map TUf,j on T*U'lj and vice veisa.-T*W^f on 
TW^j. These mappings are of the following form: 

As on a manifold of finite dimension, on there exists the unique Levi-Civita 
connection, i.e. a Riemannian one without torsion. In this case as it is shown 
in 

T[g, 0] (fc, x; OJ, v;h,e)^- I {-tr{k)tr{uj)tr{h) 

* J M 

+3tr{k)tr{LLi x h) + tr(uj)tr{k x h) + tr(h)tr{k x us) 

-Ur{k X u! X h)+ tr{h)x(p + tr{uj)xO - tr{k)(pe}N^^dfi{g) , 

and 

TZiem[g, (j>]{uj, ip; k, x, I, h, 9) 

= / {{l/i)tr[{kxOj-u!xk)x(hxl~lxh)] 
Jm 

+K^[tr{u} X l)tr{k xh)~ tr{h x ui)tr{k x I) 

+tr{oj X l)x0 - tr{u! x h)x(J + tr{k x h)(pa - tr{k x l)(pe]}N^^dfj,{g) , 

where 



— ——- -T, UJ = OJ — 



tr{uj) 



16(d- 1)' d " ' 

Consider now the dynamics with the following Hamiltonian in M-superspace 

W[g>,7r,p] = i^~^[5,(/)](7r,p;7r,p) + V[g,(j)] , 

where 

V[g,4>]^ I {-R{g) + \\\d4>f + Fmdti{g), 
Jm 

The dynamics given by this Hamiltonian 

V„u = -dV , 
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where V is a covariant derivative on the M-Superspace, corresponds to Ein- 
stein's equations together with the foUowing constraints equations, 

W[5,<A;7r,p] = , 

These equations constrain only initial conditions, but not dynamics. 

Now to investigate stability/instability of the solutions of the Einstein's equa- 
tions the following objects should be used [Tj, (in order to derive f2„) 

dV = iV({5* - c • Fg}N + Hess{N), {F' + A0} - g{d(l), dN)) 

and 

-Hv{h, x) = y[Kx)dy = 'P{h. X) + -DiK X) + ^f{h, x) , 

where 

Vih, x) = (Hh) + c[{triS xh)^ F'x}g - (ir(5) + F)h], 

~tr{Hess{(j)) x h) + F"x) ; 

V{h,x) =N^ {[^{/^h - a5h + Hess{tr{h)) - d(j) ® dtr{h)\* , 

~g{d<l>,5h+\dtr{h)) + ^x) ; 
U{h,x) ^ N{-Dg{Hess{N)).h,^g{dx,dN)) , 

c= (2(d-l))-' , 9/ = ff''%, 
T* ^T-c-tr{T)g , S ^ Ric - ^d(j) d(j) , 
Hessif) = , A/ = tr{Hess{f)) = -5'^V|a|6 , 

''(/i)ah = -Riema''b'''hcd - Rierrn," Jhcd + RiCahcb + Ricb^Ka ■ 

In further investigation [7| we will employ suggested method to the problem of 
stability of cosmological solutions. 

4 Conclusions 

Thus, we have proposed a geometrical approach to the investigation of Hamilto- 
nian systems on (Pseudo) Riemannian manifolds. We have derived the equation 
for the normal to the velocity of the flow invariant vector Eq. Q in terms of 
the new covariant derivative. We have introduced tensor il^ which characterises 
stability, instability, hyperbolicity of the integral curves of Hamiltonian systems. 
This approach is valid for any Hamiltonian system, while well known reduction 
to the geodesic flow requires some more conditions. Besides, there is no need to 
change metric and time parameter. This allows obtaining characteristic param- 
eters by means of physical values. As an important application we will consider 
in details the behaviour of the N-body gravitating systems and cosmological 
solutions of the Einstein's equations elsewhere [7]. 
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